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Abstract 

Based on quantum-kinetic equations, coupled spin-charge drift-diffusion equations are derived 
for a two-dimensional electron gas on a cylindrical surface. Besides the Rashba and Dresselhaus 
spin-orbit interaction, the elastic scattering on impurities, and a constant electric field are taken 
into account. From the solution of the drift-diffusion equations, a long-lived spin excitation is 
identified for spins coupled to the Rashba term on a cylinder with a given radius. The electric-field 
driven weakly damped spin waves are manifest in the components of the magnetization and have 
the potential for non-ballistic spin-device applications. 
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I. INTRODUCTION 



In the emerging field of spintronics, a main issue is the avoidance of spin randomization 
while stimulating controlled rotations of spins in a spin field effect transistor. Therefore, the 
recent proposal 1 for a spintronic device, which operates in the non-ballistic regime, received 
considerable interest. According to this design, spin relaxation is suppressed, when the 
Rashba and Dresselhaus spin-orbit interaction (SOI) constants are tuned by an external 
gate voltage so that their couplings become equal. In fact, the suppression is a consequence 
of an exact spin-rotation symmetry,-^ which leads to a strong anisotropy of the in-plane 
spin-dephasing time. At the presence of an in-plane electric field, the persistent spin helix 
is converted into a field-dependent internal eigenmode.- In close analogy to space-charge 
waves in crystals, these field mediated spin excitations can be probed by optical grating 
techniques.- Both spin and charge pattern, which are generated by polarized laser beams, 
provide the required wave vector for the excitation of internal eigenmodes. Suppressed spin- 
relaxation occurs not only in (001) GaAs/Al x Gai__ x As quantum wells with balanced Rashba 
and Dresselhaus SOI strengths but also in (110) quantum wells with Dresselhaus coupling. 5 

Most activities in the field of spintronics that are based on the Rashba and Dresselhaus 
SOI refer to a plane two-dimensional electron gas (2DEG), which is confined by a semicon- 
ductor quantum well. However, the diversity of spin-related phenomena markedly increases 
for different geometries of the 2DEG. In dependence on the curvature of the surface, in 
which the 2DEG resides, additional contributions to the SOI appear that may lead to new 
spin effects. Examples of current interest provide microtubes fabricated by exploiting the 
self-rolling mechanism of strained bilayers. These rolled-up structures exhibit pronounced 
optical resonances 6 arising from micron-sized cylindrical resonators or give rise to novel 
magnetoresistance oscillations, which were observed in the ballistic transport of electrons 
on cylindrical surfaces.- For non-ballistic spintronic device applications, the prediction of 
a conserved spin component, which arises when the Rashba coupling constant 71 equals 
the quantity h/2m*R (with R being the radius of the cylinder and m* the effective mass 
of the 2DEG) is most interesting.- The identification of this novel long-lived spin mode 
by fabricated curved samples seems to be feasible with the present-day technology-^ 1 ^ 1 ^ 
It is the aim of this paper to study the field dependence of the predicted spin helix by 
a systematic consideration of electric-field mediated eigenmodes of a spin-charge coupled 
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2DEG that is confined to a circular cylinder. Both Rashba and Dresselhaus SOIs as well as 
spin-independent elastic impurity scattering are taken into account. 



II. BASIC THEORY 

While in classical mechanics the restricted motion of particles on curved surfaces is un- 
ambiguously described by equations of motion, the quantum-mechanical study of curved 
systems starts from two different perspectives. In the first, widely used method, the three- 
dimensional Schrodinger equation is converted to its two-dimensional counterpart by an 
appropriate confining procedure.— This approach naturally accounts for the fact that the 
curved structures are embedded in a three-dimensional space, in which electric and mag- 
netic fields could be present. The second alternative description of the carrier dynamics on 
curved samples completely rests on a two-dimensional model.— In our study of SOI on the 
surface of a cylinder, we apply the widely accepted first approach, which was already used 
for studying spin effects on curved surfaces.— iI&iiLi& The second- quantized version of the 
Hamiltonian has the forrn^ 



oo 
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ak zS >{<p) 



in which Rashba and Dresselhaus contributions appear with the coupling constants 71 and 
72, respectively. The creation [a/ k s ((f)) and annihilation [ak zS (<p)] operators depend on the 
spin index s, the wave vector component k z along the axis of the cylinder, and the angle 
(p. The SOI terms include the Pauli matrices cr, the transverse momentum operator p v , 
and a matrix S that introduces off-diagonal elements with respect to the spin index. These 
quantities are defined by 
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The periodic boundary conditions on the cylinder surface are accounted for by a discrete 
Fourier transformation, which is applied in the form 



m=— oo m=— oo 



By this transformation, the projection of the total angular momentum on the cylinder axis 
appears and the Hamiltonian simplifies considerably. In addition to the SOI, both elastic 
scattering on impurities with the short-range coupling strength U and an external electric 
field E (applied along the cylinder axis) are taken into account. As it is assumed throughout 
the paper that the radius R of the cylinder is much larger than the lattice constant, we 
introduce the electron momentum vector 

k=(k v ,k z ,0), k v = (m+^j /R, (4) 

in order to express the Hamiltonian in a form that is very similar to the case of planar 
geometry. We obtain 

H = ^2e(k)a] u a hs + ^^(^i(fc) • tT aa >)a] u a ka > (5) 

fc . S fc S 

k,k' s k,s 



K=0 



with the parabolic dispersion relation 



h 2 k 2 h \ 

The main part of the SOI is included in the vector 

= ^0, ~{lik z - 72^), fc 9 ,(7i - 2^T^) - 72**^ • (7) 

The model description of spin-independent scattering on impurities in Eq. (jSJ) has the advan- 
tage of simplicity, permitting us an exact treatment of scattering in the Born approximation 
via the scattering time r defined by 

1 2irU 2 



t h 



J25(s(k)-s(k')). (8) 



fc' 

All information about the spin-orbit coupled electron ensemble is provided by the spin- 
density matrix 

f s Ak,k'\t) = (ala k , s ,) t , (9) 



which is calculated from quantum-kinetic equations. A transparent physical interpretation 
of final results is facilitated by considering the projected spin vector on a local trihedron. 
This transformation is achieved by 

f = Y,fs'Sss', (10) 

s,s' 

with the following transformation matrices 

< sv =li ° ie2liP I - s ' ; = o I I- < s "' -. \ , I an 





2 I -ie" 2 ^ 

which project to the cylinder axis (S z ), as well as to the tangential (S^) and normal (S r ) 
directions. To proceed, the wave vectors are shifted according to k — > k + re/2 and k' — > 
k — n/2 with k v = (m + m' + 1)/2R and = (m — m') / R. The derivation of spin-charge 
coupled kinetic equations is carried out by applying the same steps as in our previous study 
of the planar geometry— The final result is expressed by the kinetic equations 

-f(k, n\t) -—k-Kf + iu{n)-f+ e -E- V fc / = -(7 - /), (12) 
at m* n t 

-f(k, K \t)-—(k- K)f - 2oj(k) xf + + Ue ■ V fc )/ 

ot m* n 

in which a new SOI vector appears given by 

iv (k) = sin(2</?), uJi v (k) cos(2v9), -^(k)). (14) 

On the right-hand side of Eq. ( |T3i) . there remain scattering contributions that are propor- 
tional to the SOI. These terms are necessary for a consistent treatment of a homogeneous 
2DEG. The bar over quantities in Eqs. (|T2l) and (|T3l) indicates an integration over the polar 
angle a of the vector k = fc(cos a, sin a, 0). The kinetic Eqs. f|T2|) and f|T3|) serve as a starting 
point for various studies of spin effects on cylindrical surfaces. We mention only the field- 
induced spin accumulation and the influence of the spin degree of freedom on the charge 
current (as well as the appearance of a "spin current" on the cylinder).— In this paper, we 
do not follow this interesting line of reasoning but look for a solution of Eqs. (fT2|) and (|T3|) 
in the long-wavelength and low-frequency drift-diffusion regime. 



The envisaged macroscopic behavior of the coupled spin-charge system is established 
during an evolution period, in which a nonequilibrium spin polarization and charge density 
still exist, whereas the energy of particles already thermalized.— In this transport regime, 
the following separation ansatz for the mean components / and / is justified 

n'(e(k)) 



f(k, n\t) = -F( K \t)^4^-, f(k, K \t) = -F(K\t)- 



(15) 



dn/deF dn/deF 
where n(e(k)) denotes the Fermi distribution function and n = J dep(e)n(e) is the equilib- 
rium carrier density (with p(e) being the density of states). n'(e(k)) is a short-hand notation 
for dn/deik). Adopting this approximation also for the field contributions on the left-hand 
side of Eqs. ( I12p and fll3p . we obtain a set of linear equations for the components of the 
spin-density matrix. The solution is expanded with respect to n and integrated over the 
angle a. This procedure leads to coupled equations for the charge F and spin F distribu- 
tion functions, the solution of which is easily integrated over the energy e[k). The resulting 
coupled spin-charge drift- diffusion equations take the form 



d_ 

dt 



- ijxE ■ k + Dk 2 



7 Z/iTD T 

F + — uj {k)-M + — — ([Ax«(k)].M) = 0, (16) 



d_ 

dt 



ipE ■ n + Dk 2 + T 



M H —M x H eff 



eff) 



-2i — -[A x uj(k)]F — G, 



(17) 



n c 

for the charge density F and the magnetization M = fisF with = eh/2m*c being the 
Bohr magneton. The vector G on the right-hand side of Eq. (fl7|) accounts for the source 
of an external spin generation. Furthermore, D and \i denote the diffusion coefficient and 
the mobility that are related to each other via the Einstein relation /x = eDn'/ n. The Pauli 
susceptibility is given by x — ■ Scattering times that refer to various spin components 
are collected by the symmetric matrix T, which is given by 

ADm* 2 



h 2 



+ aj 2 + alx + a| 2 — (022032 + «2ia3i) — (ana2i + 022^12) 
(022032 + 021031) a\i + o? 2 + a li + a 22 -(012O32 + a n a 3 i) 



\ 



a 21 ~l~ a 22 "I" O31 



— (ana2i + 022O12) — (012O32 + 011031 
where the quantities expressed by the spin-orbit coupling constants 



' 32 



(18) 



/ 



on = 7 2 sin(2vj), a 2 i = 7 2 cos(2y2), a 12 = -7isin(2<^), a 22 = -7iCOs(2y2), (19) 



The electric field is accounted for by the vector 



Q32 = 72- 



(20) 



A = (a 2 i/jE v + a 22 fxE z , a^iiE^ + a 32 fxE z , a\\[iE v + ai 2 fxE z ), (21) 

from which an effective magnetic field 

H eff = (A + 2iDlj(k)), (22) 

is derived that enters Eq. (fTTj) for the field-induced magnetization. The appearance of a 
magnetic field H e ff, which is solely due to the electric field, illustrates why there is a 
perfect electric-field analog of the Hanle effect.— 

The drift-diffusion Eqs. (1161) and (fT7|) for the charge density F(k, t) and magnetization 
F(K,,t) provide the basis for the study of many spin-related phenomena of a curved 2DEG 
in the drift-diffusion regime. Here, we shall focus on spin-related eigenmodes of the system. 



III. LONG-LIVED SPIN WAVES 

A solution of Eq. ffTTI) is searched for under the condition that the retroaction of spin on 
the charge density can be neglected so that the carrier density is given by its equilibrium 
value (F = n). Performing a Laplace transformation with respect to the time variable t and 
introducing the abbreviations 

M' = M- xH eff , E = s - ifiE ■ k + Dk 2 , (23) 

with s being the Laplace variable, Eq. fflTj) is converted into the linear equations, 



EM' + YM' + —M' x H eff (24) 

m*c 

= M(0) + G/s - 2i — "A x uj(K)n - x^H eff , 

which are symbolically written as TM' = Q. Eigenmodes of the spin subsystem are calcu- 
lated from the zeros of the determinant of the matrix T. A simple but cumbersome algebra 
leads to the result 

detf = E(a 2 + «|) + g 2 (a + ^) , (25) 
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in which the short-hand notations ojh = { e / m * c )H e ff and o = S+gi are used. The coupling 
constants g\ and g 2 are given by 
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7i +7 2 
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2m*R 
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7i 
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92 



ADm 



*2\ 2 



7 2 - 7i 7i - 



h 



2m*R 



(26) 



(27) 



The cubic equation det T = with respect to the Laplace variable s has three solutions, 
which give the dispersion relations of spin excitations. Most eigenmodes have a finite lifetime. 
However, there is one long-lived spin excitation, whose damping completely disappears for 
a given wave number k z . This mode appears for a model without any Dresselhaus SOI 
(72 = 0), when the coupling constant 71 matches the quantity h/2m*R. In this case, we 
obtain (s — > iuj) 

uj 1>2 = -nE z (K z ±K)-iD(K z ±K) 2 , (28) 

with K = 2m*^fi/h being a wave number that is built from the Rashba spin-orbit coupling 
constant 71. This soft mode becomes increasingly undamped in the limit k z —>■ K. The 
persistent spin mode of this kind, which is a consequence of a new spin-rotation symmetry, 
has no counterpart in the planar Rashba model and is a distinct feature that solely appears 
on a cylinder surface. 

In order to excite the persistent spin wave, a regular lattice of spin polarization Q r 
perpendicular to the cylinder surface is provided by laser pulses. For simplicity, the spin 
generation is assumed to have the form 

QrO 1 



Qr 



■[8(k z - Kq) + 5(k z + Kq)]. 



(29) 



Under the condition Q v = Q z = 0, the solution M = xH e ff + T l Q of Eq. (fT7l) is expressed 
by 



M z 



liE z K 



9 1 2 v ^ r ' 



Q r , M r 



a 



-Qr 



(30) 

In the derivation of these equations, it was considered that the inverse Fourier transformation 
with respect to k v leads to (p — 0. The inverse Laplace transformation and the integration 
over k z give for the non- vanishing components of the field-mediated magnetization the final 
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results 



M r (z,t) 



QrO 



e -D(K +K) 2 t 



COs[k,qZ + fiE z (ri Q + K)t] 



2 



+e -D(K -K)H cos { KqZ + j^E^Ko - K)t] 
M z (z,t) = Mi-\z,t)-Mi + \z,t), 



(32) 



(31) 



with 



MT(z,t) 




2Dkq cos[k z + [iE z (k ± K)t] 



D(K ±K) 2 t 



(33) 



Both components M z and M r consist of a strongly and weakly damped oscillating term. 
Under the resonance condition k = K, the first mode quickly disappears, whereas the 
second mode becomes completely undamped. A smooth dependence on the electric field E z 
persists in the magnetization M z along the cylinder axis. A slight detuning of the resonance, 
however, leads to the appearance of an electric-field driven spin wave, the damping of which 
is extremely weak. The frequency of this long-lived spin excitation is directly controlled by 
the applied electric field. The situation is similar to the persistent spin helix of a planar 
2DEG. Therefore, it is supposed that the robust spin wave on a cylinder and its direct 
manipulation by an electric field has the potential to be utilized in future spintronic device 
applications. 

IV. SUMMARY 

Nanostructures with a great variety of novel geometries like curved graphene systems and 
rolled-up 2DEG are now experimentally available. Hence, the rigorous theoretical study 
of the dynamics on such curved surfaces became a subject of recent interest. Especially 
spin effects have been treated because the curvature of the surface gives rise to additional 
contributions to the SOI. Consequently, the number of possible spin effects considerably 
increases in nanostructures with curved geometries. This observation further stimulates 
activities in the field of spintronics. A key point regarding spin-field-effect transistors is 
the exclusive manipulation of spin by means of an electric field. Particularly attractive is 
the proposal for a device working in the non-ballistic regime, where spin scattering in a 
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planar 2DEG is suppressed due to a spin-rotation symmetry.- A similar effect that occurs 
on the surface of a cylinder was studied in the present paper. Based on quantum-kinetic 
equations for the spin-density matrix, rigorous coupled spin-charge drift-diffusion equations 
were systematically derived for a cylinder, whose radius is much larger than the lattice 
constant. From the solution of these equations, the dispersion relation of field-dependent 
spin eigenmodes are identified. In general, there are three damped spin excitations, the 
character of which is determined by the coupling constants 71 and 72 of the Rashba and 
Dresselhaus SOI. For the pure Rashba model (72 = 0), a long-lived spin wave exists, when 
the radius R of the cylinder matches the condition R = H/2m* r y\. This finding is of particular 
interest as an applied electric field stimulates a nearly undamped spin wave. The excitation 
mechanism of spin waves has the same character as the excitation of space-charge waves, 
which are normally strongly damped. Nevertheless, space-charge waves in crystals were 
clearly demonstrated in experiment. To my knowledge, completely undamped space-charge 
waves do not exist. Their damping is reduced, however, in the regime of negative differential 
conductivity due to a negative Maxwellian relaxation time.— The complete disappearance 
of the damping of an excitation is a novelty that occurs in special spin subsystems with a 
fc-linear SOI. The above mentioned peculiarity of the Rashba model on a cylindrical surface 
has no counterpart in a planar 2DEG. Unfortunately, the experimental demonstration of 
this effect is rendered more difficult because the huge internal strain within the tube breaks 
the bulk inversion symmetry so that an appreciable Dresselhaus contribution to the SOI is 
expected, which detunes the strong spin resonance. If this problem can be circumvented, 
the long-lived field-mediated spin excitations on a cylinder have the potential to be utilized 
in spintronic devices that work even in the non-ballistic regime. 
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